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, ronments over time. By means of a multi-fluid model that treats plasma and neutral hydrogen 



At present, the heliosphere is embedded in a warm low density interstellar cloud that belongs 
to a cloud system flowing through the local standard of rest with a velocity near ~18 km s^^. The 
velocity structure of the nearest interstellar material (ISM), combined with theoretical models 
of the local interstellar cloud (LIC), suggest that the Sun passes through cloudlets on timescales 
of < 10^-10'* yr, so the heliosphere has been, and will be, exposed to different interstellar envi- 



self-consistently, the interaction of the solar wind with a variety of partially ionized ISM is in- 
vestigated, with the focus on low density cloudlets such as are currently near the Sun. Under 
the assumption that the basic solar wind parameters remain/were as they are today, a range of 
ISM parameters (from cold neutral to hot ionized, with various densities and velocities) is consid- 
ered. In response to different interstellar boundary conditions, the heliospheric size and structure 
change, as does the abundance of interstellar and secondary neutrals in the inner heliosphere, 
and the cosmic ray level in the vicinity of Earth. Some empirical relations between interstellar 
parameters and heliospheric boundary locations, as well as neutral densities, are extracted from 
the models. 

Subject headings: cosmic rays — hydrodynamics — interplanetary medium — ISM: clouds — ISM: 
structure — stars: winds, outflows 

1. INTRODUCTION 

by the solar wind - LISM interaction, and they 
The heliosphere is a low density cavity that vary in response to the Galactic environment of 
is carved out from the local interstellar medium the Sun as the Sun and interstellar clouds move 
(LISM) by the solar wind. The size and par- through space. The path of the Sun has taken us 
tide content of the heliosphere are determined through the Local Bubble void (galactic longitudes 
180° < I < 270°, Frisch & York 1986), and we 

^also at: Institute of Geophysics and Planetary Physics, University of CaUfornia, Riverside, CA 92521. 
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have recently ( ^ 10^ — lO'"^ yr ago, depending on 
cloud shapes and densities) entered a clumpy flow 
of low density interstellar material (Prisch 1994). 
This clumpy flow, the "cluster of local interstel- 
lar cloudlets" (CLIC), is flowing away from the 
Sco-Cen association and extends 10-30 pc into the 
Galactic center hemisphere and ^ 3 pc for many 
directions in the anticenter hemisphere. Inhomo- 
geneities in the CLIC create temporal variations in 
the dynamic interstellar pressure at the solar lo- 
cation, which may produce significant variations 
in heliosphere properties over geologically short 
timescales (Frisch 1993, 1997; Zank & Frisch 1999; 
Florinski et al. 2003a; Frisch et al. 2005; Frisch 
2004). 

The heliosphere itself is a dynamically chang- 
ing object which is highly sensitive to interstel- 
lar pressure (e.g. Holzer 1989; Zank 1999). The 
interaction of the ISM with the fully ionized so- 
lar wind gives rise to the heliospheric morphology 
which includes the heliopause (HP), a tangential 
discontinuity separating solar wind and LISM, and 
the termination shock (TS) where the solar wind 
becomes subsonic and is diverted downstream to 
form a heliotail. Depending on the pressure of 
the surrounding interstellar material, an interstel- 
lar bow shock (BS) may form upwind of the he- 
liopause. These general boundaries are created by 
the plasma interaction, yet the presence of neu- 
tral H and its coupling to the plasma protons 
via charge exchange greatly influences the details 
of the heliospheric morphology and the location 
of its boundaries (see Zank (1999) for a review). 
The sensitivity of the heliosphere to variations in 
the physical characteristics of the interstellar cloud 
surrounding the solar system is poorly understood, 
and in this paper we focus on heliosphere varia- 
tions due to encounters with a range of low den- 
sity clouds such as expected in the immediate past 
and future solar history. 

Different interstellar environments may pro- 
duce noticeable changes in the interplanetary envi- 
ronment of the inner heliosphere, as indicated by 
the amount of neutral H, anomalous, and galac- 
tic cosmic rays (GCR) at 1 AU. There is some 
evidence that lunar soils contain an archive of iso- 
topic abundances that are different from the par- 
ticle environment of the present era (Wimmer- 
Schweingruber & Bochsler 2000), and antarctic ice 
cores show signatures that may be interpreted as 



cosmic ray background variations at Earth (Rais- 
beck et al. 1987; Sonett et al. 1987; Florinski & 
Zank 2006). These possibilities have motivated 
our study of the behavior of the global heliosphere 
under variable boundary conditions resulting from 
passage through interstellar clouds. 

Given the inhomogeneity of the local solar 
neighborhood and the galactic environment in 
general, we test the heliosphere response to a 
range of local interstellar boundary parameters 
using about two dozen specific parameter sets. 
Our choices are justified in §2. Four highlights 
of the corresponding heliospheric models are de- 
tailed in §3. The results of all the heliospheric 
models calculated for this study suggest relation- 
ships of the heliospheric boundary locations and 
the neutral particle densities with the interstellar 
parameters, discussed in §4. The synopsis of all 
the individual model results through these rela- 
tions is a quantitative expression of the sensitivity 
of the heliosphere to changing interstellar bound- 
ary conditions. Furthermore, the relations allow 
for a prediction of boundary locations and par- 
ticle content for heliospheres with yet different 
boundary parameters, without actually engaging 
in a complex, non-linear global heliosphere simu- 
lation. This predictive power can also be used in 
the emerging field of astrospheres, which are the 
analogues of heliospheres around solar-like cool 
stars. 

We discuss the response of the global helio- 
sphere to variable interstellar properties, and spec- 
ulate on aspects of the implications of these vari- 
ations for the 1 AU location of the Earth in §5. 

2. PROPERTIES OF ISM IN THE SO- 
LAR NEIGHBORHOOD 

The Sun is embedded in a flow of warm low 

density gas with an upwind direction in the local 
standard of rest directed towards the Scorpius- 
Centaurus Association (Prisch 2004). The CLIC 
is defined by high resolution absorption lines to- 
wards nearby stars, and lower resolution observa- 
tions of white dwarf stars in the far ultraviolet 
(UV) and extreme UV. With the possible excep- 
tion of a Oph, interstellar column densities to- 
wards stars within 35 pc of the Sun do not ex- 
ceed cm-2 (e.g. Frisch 2004; Redfield & 
Linsky 2004a; Wood et al. 2005b). To date, over 
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150 absorption components have been identified 
by velocity in optical and UV observations of at 
least 90 stars sampling the nearby ISM (see refer- 
ences in Frisch ct al. 2002; Rcdficld & Linsky 2002, 
2004a,b). The ISM towards nearby white dwarf 
stars is partially ionized (Vallerga 1998; Holberg 
et al. 1999; Frisch 2004), and local variations in 
iV(Fe+)/iV(D°) show that the CLIC is inhomoge- 
neous (Frisch 2004). Prom these data, a general 
picture has emerged that the ISM within ^35 pc 
of the Sun is dominated by low density warm gas. 
The general properties of this nearby ISM are con- 
sistent with partially ionized, low column density 
gas (log A^(H°) < 18 dex, N in cm~^) described by 
radiative transfer models (Slavin & Frisch 2002). 

2.1. Short-term Variations in the Solar 
Environment 

The properties of the CLIC are diagnosed by 

Dopplcr-broadencd absorption features represent- 
ing clouds (or "cloudlets" ) observed in the optical 
and UV data towards ~100 nearby stars. The 
best-fitting flow velocity in the Local Standard of 
Rest (LSR) is — 19.4±4.6 km s~^, with an upwind 
direction I = 331.4°, b = -4.9°. For comparison, 
the LSR local interstellar cloud (LIC) velocity is 
-20.6 km s~^, and the upwind direction is {I, b) 
= (317.8°, -0.5°). These LSR values assume the 
standard solar apex motion of 19.7 km towards 
I = 57°, b ~ +22°. The corresponding heliocentric 
flow vector is 28.1±4.6 km s~^ from the upwind 
direction {I, b) = (12.4°, 11. 6°). ^ The flow velocity 
is somewhat sensitive to the star sample because a 
velocity gradient between the upwind and down- 
wind directions indicates the flow is decelerating 
(Frisch & Slavin 2006). 

Data for the nearest cloudlets have been pre- 
sented in a scries of studies by Rcdficld & Linsky 
(2002, 2004a,b). Focusing only on UV observa- 
tions of cloudlets within 15 pc as a predictor of 
past and future variations in the Galactic envi- 
ronment, we find a range of temperatures T and 
turbulent velocities ^,T= 1700-12,600 K and ^ = 
0-5.5 km s~^, with a mean temperature 6780±190 



^An alternate solar apex motion, based on Hipparcos data 
(13.4 km s-i towards I = 27.7°, b = 32.4°, Dehnen & 
Binney 1998), yields an LSR bulk flow velocity -17.0±4.6 
km with upwind direction {I, b) = (2.3°, -5.2°) and an 
LSR LIC vector of -15.7 km s-\ upwind {I, b) = (346.0°, 
0.1°). 



K, comparing favorably to the LIC temperature 
6300±340 K inferred by spacecraft (Witte 2004). 
Only ~35% of space within 10 pc of the Sun is 
filled with neutral gas if this material has the same 
density as the LIC (~0.20 cm~^, Frisch & Slavin 
2003), and if iV(D°)/iV(H°) = 1.5 x 10"^ The 
mean cloud lengths arc 0.9±0.3 pc. At a relative 
Sun-cloud velocity of 19 km s~^, such a distance 
is traversed in ~47,000 yr. The ISM filling factor 
/ found locally varies from / ^ 0.60 towards a 
Aql (/, b, d = 47.7°, -8.9°, 5.1 pc), to / ~ 0.29 
in the opposite direction towards Sirius {I, b, d = 
227.2°, -8.9°, 2.6 pc). Here, / is the fraction of 
space filled with ISM if all ISM has a density of 
0.2 cm"^. 

Several sets of data indicate that the ISM 
within ^3-10 pc is not uniform. The ratio 
A^(Fe"'")/iV(D°) varies between the downwind and 
upwind LSR directions by up to a factor of ^8, 
apparently from ionization or abundance varia- 
tions (Frisch 2004). The ISM temperature within 
5 pc varies by over a factor of 4 (Redfield & Linsky 
2004b). If the cloud in front of the nearest star a 
Cen also extends in front of a Oph, as indicated 
by their common velocity, and is uniform, then 
Ca+ and 11° data suggest a density n > 5 cm~^ 
(Frisch 2003), in contrast to the LIC density ~0.3 
cm~^. Velocity variations of 10 km s~^ or more 
are also found along several individual sightlines 
(Frisch et al. 2002; Redfield & Linsky 2004b). 

A series of radiative transfer models appropri- 
ate for the radiation field and physical properties 
of the low column density material close to the Sun 
show that equilibrium occurs for a range of ioniza- 
tion levels in low density ISM (Slavin & Frisch 
2002; Frisch & Slavin 2006). Furthermore, the 
boundary conditions of the hcliosphcre will vary 
as it traverses low density ISM strictly because of 
ionization variations within the cloud. These mod- 
els generally consider clouds with ritot ^ 0.3 cm~^, 
iV(H°) < 10^* cm^^, a local radiation field consis- 
tent with observations of the diffuse radiation field 
between 3000 A and 0.5 keV at the solar location, 
and allow for additional radiation emitted by a 
possibly magnetized conductive interface between 
the local warm gas and adjacent very hot Local 
Bubble plasma. These models predict ISM equi- 
librium conditions for n(II°) = 0.16-0.26 cm~'^, 
hydrogen ionization levels x(H) = II+/(H°+II+) 
= 0.19-0.34, and cloud temperatures of 4900-8300 
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K. 

The time of the Sun's entry into, and exit from, 
the Lie can be estimated using observations of H° 
and D° towards nearby stars combined with Ti(H°) 
derived from the Slavin & Frisch radiative transfer 
models. These models of the LIC indicate that 
n(H°) = 0.19 0.21 cm-3 and n(H+) - 0.1 cm-^ 
at the solar location, and that n(H°) decreases by 
<20% between the Sun and surface of the LIC. 
Assuming a constant LIC density of n(H°) = 0.2 
cm~'^ and using the limits on the LIC component 
towards 36 Oph (A^(H°) < 6 x 10^^ cm'^, Wood 
et al. 2000a), we infer that the distance to the LIC 
surface in this direction, as defined by a velocity 
discontinuity in the gas, is <0.1 pc, suggesting the 
Sun will exit the LIC in less than 3700 yrs. 

The entry of the Sun into the LIC can be calcu- 
lated after transforming into the LSR frame and 
assuming a LIC morphology (e.g. Prisch 1994). If 
the LIC velocity vector is perpendicular to the 
surface, and assuming A'^(H°) = 4.0 x 10^^ cm~^ 
towards a CMa (Hebrard et al. 1999), the Sun 
will have entered the LIC ^6700/11,500 yr ago 
for the Hipparcos/Standard solar apex motion, re- 
spectively. The assumed column density requires 
additional ISM near the LIC velocity towards the 
downwind stars a Aur and Ori. Alterna- 
tively, the LIC column densities towards down- 
wind stars can be used to define a plane that ad- 
vances through space with the LIC velocity vec- 
tor. The useful downwind stars for this estimate 
are a CMa, a CMi, and a Aur, where 

column densities for the LIC are, respectively, log 
iV(H°)= 17.60, 17.90, 17.80, 18.26 cm^^ based on 
D/H = 1.5 X 10^'' and data by Hebrard et al. 
(1999) and Rcdficld & Linsky (2004a); the selec- 
tion of any three of these four stars, and assuming 
n(H°) = 0.2 cm^'^, suggests that the Sun entered 
the LIC about 28,000-30,250 yr ago. When un- 
certainties of ~30% are incorporated to reflect the 
various assumptions, these estimates suggest that 
the Sun has entered the interstellar cloud compo- 
nent at the LIC velocity sometime within the past 
40,000 yr, and will exit it sometime within the 
next 4000 yr. 

The star Ori is within 15° of the downwind 
direction, and shows a cloud with a relative Sun- 
cloud velocity of 21.6 km s^^. The Sun and this 
cloud would have first crossed paths -^47,000 yr 
ago for n(H°) = 0.2 cm~^. Beyond x^ Ori) the 



next neutral gas in the downwind direction is over 
50 pc away. Allowing for ~30% uncertainties and 
possible gaps between clouds, the Sun would have 
entered the CLIC within the past ~60,000// yr. 

2.2. Variations in the Global ISM 

An ISM with a wide range of properties is found 

within ~350 pc of the Sun. including low density 
hot gas in the Local Bubble that emits soft X-rays. 
UV and radio observations of low column density 
ISM show that a range of ISM types arc possible 
at low column densities. The ISM within that dis- 
tance provides a model for the types of ISM the 
Sun may encounter over timoscalcs of several mil- 
lion years. The Sun will move lij 16-20 pc through 
the LSR per million years, and interstellar clouds 
(with velocities of up to ^ 100 km s~^) may move 
hundreds of parsecs.^ 

One example for ISM structure in this range is 
provided by Welty ct al. (1999) who compare opti- 
cal and UV ISM data for the star 23 Ori (300 pc), 
which is in the direction of the Orion-Eridanus 
soft X-ray supcrbubblc, and find a complex sys- 
tem of cloudlets showing a wide range of proper- 
ties, representing typical diffuse ISM. Twenty-one 
cloudlets with LSR velocities that range from 120 
to -|-8 km s~^ are found. Four low velocity clouds 
(positive Vlsr) at ~100 K are present. They 
are massive (log A^(II) = 20.7 cm"2)^ moderately 
dense (10-15 cm~'^), suggesting a cloud thickness 
of ^ 15 pc, and primarily neutral. Crossing such 
a cloud might take the Sim ^1 Myrs. Warmer 
dense clouds (~ 15-20 cm~^, primarily neutral, 
~3000 K) with thicknesses ^ 1 pc are found at 
moderately low velocities (Vlsr ~ — 17 ^ km 
s~^). At higher velocities (Vlsr ~ —60 -17 km 
s~^) low density gas (n(H°) = 1-6 cm~^) is found. 
This material appears to be in thin sheets with 
thicknesses of 0.001-0.04 pc. Warm (~8000±2000 
K), rapidly moving (Vlsr = —130 — > -100 km 
s""'^) shocked low column density clouds are also 
seen. This gas is partially ionized (n(e~) = n(H°) 
= 0.4-0.5 cm~^) and arises in clouds with thick- 
nesses of about 0.005 0.12 pc. The ionization of 
this high velocity gas suggests an interstellar ra- 
diative shock where the gas is not in ionization 
equilibrium (since the coUisionally ionized species 
show T ~ 25,000 K, but Doppler ^-values indicate 



velocity of 1 km s ^ corresponds to ~1 pc/Myrs. 
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Fig. 1. — Velocity distribution of the CLIC versus 

CNM and WNM: The CLIC components are plot- 
ted with the symbol "x". The CLIC components 
shown here are towards stars within 50 pc. The 
LSR velocities arc derived from 11°, D°, and op- 
tical Ca+ data in Frisch et al. (2002); Redfield & 
Linsky (2002, 2004a); Wood et al. (2005b) (Stan- 
dard solar apex motion used). Points based on 
iV(D°) (or iV(Ca+)) data assume A^(D°)/iV(H°)= 
1.5x10-5 (or A/'(Ca+)/Ar(H°)= 1.0x10-^). Filled 
and open circles show the CNM and WNM 11° 
components, respectively (from Heiles & Troland 
2003a). Only high latitude HI sightlines are in- 
cluded (|b| >25°). Note that the CNM and WNM 
data sample may contain non-local components 
with velocities affected by differential galactic ro- 
tation. By selecting the high latitude HI sight- 
lines, the sample is more likely restricted to nearby 
regions, where galactic rotation effects make mini- 
mal contributions to the CNM and WNM compo- 
nent velocities. 



6000-12,000 K). 

The Millennium Arecibo 11° 21-cm radio survey 
of warm neutral material (WNM) and cold neu- 
tral material (CNM) also provides a comparison 
sample for the CLIC (Heiles & Troland 2003a,b). 
Figure 1 shows the LSR velocities of CLIC com- 
ponents observed in the optical and UV (data 
from Frisch et al. 2002; Redfield & Linsky 2004a,b; 
Wood et al. 2005b) compared to CNM and WNM 
velocities. Except for extra high-latitude infalling 
ISM flows at V < —25 km s"^ (Lockman & 
Gehman 1991), the kinematics of the CLIC are 
similar to WNM and CNM. If viewed from the 
outside, the CLIC would appear as a medium ve- 
locity flow (17-20 km s"^) with low column densi- 
ties (iV(H°) < 10^9 cm-2). The WNM has upper- 
limit kinetic temperatures of 500 K to over 10,000 
K, and median column densities of l.Sx 10^" cm~^. 
The dominance of low mass warm clouds at inter- 
mediate velocities (<-17 km s~^) seen in Fig. 1, 
combined with the fact that ~60% of the H° is 
WNM, suggest that warm low density clouds are 
the most likely to be encountered by the Sun over 
the next million years. 

The recent discovery of cold (<100 K) tiny neu- 
tral cloudlets in the ISM, iV(H°) ~ lO^^ cm'^, 
including one towards 3C286 at a velocity within 
1.5 km s~^ of the G-cloud velocity (Stanimirovic 
& Heiles 2005), shows that tiny cold neutral clouds 
are widespread but infrequent. CNM components 
with A^(H°) < 10^* cm-2 g^j^^j densities ~20 cm'^, 
similar to values found by Welty et al. towards 23 
Ori, would have thicknesses of <0.02 pc and if 
at rest in the LSR would perturb the heliosphere 
boundary conditions on timescales of ^^100 yr. 

2.3. Modeled Clouds 

The above discussion of the ISM in the solar 
neighborhood, and close to the Sun, provides the 
basis for selecting a representative set of boundary 
conditions for modeling. The heliosphere configu- 
ration has been modeled for 27 cloud types with 
densities varying from 0.005-15 cm~^, ionizations 
ranging up to 100%, and relative Sun-cloud veloc- 
ities of up to 100 km s~^. The boundary param- 
eters are listed in Table 1. Most of the assumed 
cloud types are warm and low density clouds, but 
the possible velocities vary by an order of magni- 
tude. The Sun moves through the LSR at ~14- 
19.5 km s~^, while warm diffuse clouds have LSR 
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Table 1 
Model Boundary Parameters 



# 




nH+ 


'^tot 


X(H) 


t'LISM 


Ilism 




(cm ") 


(cm ") 


(cm '*) 




(km s ) 


(K) 


1 


0.00 


0.005 


0.005 


1.00 


13.4 


1260000 


2 


0.14 


0.10 


0.24 


0.42 


8.3 


7000 


3 


0.24 


0.04 


0.28 


0.14 


15 


3000 


4 


0.24 


0.04 


0.28 


0.14 


26 


7000 


5 


0.216 


0.047 


0.26 


0.18 


26 


7000 


6 


0.242 


0.074 


0.32 


0.23 


26 


7000 


7 


0.24 


0.10 


0.34 


0.29 


26 


7000 


8 


0.235 


0.106 


0.34 


0.31 


26 


7000 


9 


0.14 


0.10 


0.24 


0.42 


25 


5650 


10 


0.14 


0.10 


0.24 


0.42 


26 


8000 


11 


0.04 


0.04 


0.08 


0.50 


15 


3000 


12 


0.04 


0.04 


0.08 


0.50 


26 


7000 


13 


0.10 


0.10 


0.20 


0.50 


26 


8000 


14 


0.40 


0.40 


0.80 


0.50 


100 


8000 


15 


11.00 


0.15 


11.15 


0.01 


26 


100 


16 


15.00 


0.20 


15.20 


0.01 


26 


10 


17 


15.00 


0.20 


15.20 


0.01 


26 


3000 


18 


0.14 


0.10 


0.24 


0.42 


31.5 


5650 


19 


0.14 


0.10 


0.24 


0.42 


37.7 


5650 


20 


0.14 


0.10 


0.24 


0.42 


45.2 


7000 


21 


0.14 


0.10 


0.24 


0.42 


50.8 


7000 


22*^ 


0.14 


0.10 


0.24 


0.42 


68 


8000 


25 


0.14 


0.10 


0.24 


0.42 


86 


8000 


26 


0.24 


0.04 


0.28 


0.14 


50.8 


10 


27 


0.96 


0.04 


1.00 


0.04 


50.8 


10 



"^Model 23, Tsw= 2 x 10^ K. Model 24, t;sw= 500 km 3"^ 
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velocities 20 60 km s~^, and higher. In particular, 
the cloud cluster near the Sun shows depletions 
characteristic of shocked interstellar gas (Frisch 
1981, 2004), where large peculiar motions might 
be expected. Hence relative Sun-cloud velocities 
may range over an order of magnitude, leading to 
appreciable variations in the heliosphere morphol- 
ogy- 
Table 1 gives the LISM boundary conditions 
as neutral hydrogen number density n(H°), pro- 
ton number density n(H"'"), and heliocentric ve- 
locity i^LiSM and temperature Tlism of the inter- 
stellar wind, ntot is the total hydrogen density, 
and x(H) = n(H+)/ntot the interstellar hydrogen 
ionization fraction. Model 1 (ntot=0.005 cm~^, 
X(H)=1, u=13.4 km s^^, log T=6.1 K) represents 
the hot plasma interior of the Local Bubble (Snow- 
den et al. 1997). Model 2 also tests subsonic in- 
terstellar conditions, but with a 42% partially ion- 
ized dense, warm ISM (ntot=0.24 cm-^, T=7000 
K). The velocity in model 2 (8.3 km s~^) is com- 
parable to the heliocentric velocity of the blue- 
shifted cloud found towards e CMa and a CMa 
(Gry & Jenkins 2001), which may have been the 
first warm cloud encountered by the Sun as it ex- 
ited the plasma interior of the Local Bubble. 

Models 3-10 test the effect of small density 
('^tot =0.24-0.34 cm~^), warm temperature (3000- 
8000 K), and velocity (15-26 km s~^) variations 
on the heliosphere morphology, with varying ion- 
izations (x(II)=0.14-0.42), representing the warm 
low density ISM described earlier. Neutral clouds 
with T ~ 3000 K are widespread and evidently 
thermally unstable in the absence of magnetic 
pressure (Hcilcs 2001). Models 4-10, with v - 26 
km s^^, represent variants of the contemporary 
heliosphere whose interstellar boundary densities 
fall within the constraints of the observations. 
Model 9 corresponds to the a Cen environment 
(Linsky & Wood 1996; Wood et al. 2001), pro- 
vided the ISM spreads uniformly throughout the 
sightline towards this nearest star. Model 10 is 
based on the local cloud towards the white dwarf 
star RE J 1032+532 (Holberg et al. 1999). 

Models 11-14 represent the effect on the helio- 
sphere of density (ntot=0.08-0.8 cm"^) and veloc- 
ity (15-100 km s"-'^) variations in 50% ionized in- 
terstellar hydrogen at warm temperatures (3000- 
8000 K), such as might be expected for kinemati- 
cally perturbed low density gas. Models 15-17 test 



the expected dramatic differences in heliosphere 
configuration anticipated from the encounter with 
a denser (ntot=ll-15 cm~^) neutral (x(H)=0.01) 
cold or tepid (T=10 3000 K) interstellar cloud at 
the Lie velocity. Model 15 is based on the strong 
line, low velocity gas towards 23 Ori (Welty et al. 
1999), while model 17 is based on both the warm 
low velocity gas towards 23 Ori, and the thermally 
unstable warm H° gas observed at 21 cm (Heiles 
2001). Model 16 was calculated solely for compar- 
ison with model 17, representing a temperature re- 
duction by a factor of 300 without a correspond- 
ing density increase that typical Galactic values 
of the roughly constant product ntot Tlism would 
suggest. 

Models 18-25, together with 2, test the param- 
eter space around a partially ionized, diffuse inter- 
stellar cloud (x(H)=0.42; ntot=0.24 cm^-''), sam- 
pling velocity variations from 8 to 100 km s~^. 
Changes that arise from small variations in the 
solar wind are also considered (models 23 and 24). 
Model 14 is an example of cooled high velocity 
shocked gas such as a superbubble shell formed 
from a supernova shock sweeping up the ISM in 
a starforming region. Its parameters are based 
on high velocity gas observed towards 23 Ori and 
C Ori (Welty et al. 1999, 2002). Finally, models 
26 27 test a cold (10 K), high velocity (51 km 
s~^) regime, again at lower values of ntot Tlism 
than typically encountered, in order to discuss 
low-temperature regimes without getting into the 
complications of very high density heliospheres. 

3. INDIVIDUAL MODEL RESULTS 

To characterize the large-scale heliospheres that 
result when the Sun is embedded in the different 
parts of the ISM as described above, we make 
use of the multi-fluid model developed by Zank 
et al. (1996b). The multi- fluid code simultane- 
ously solves the time evolution of four interpene- 
trating fluids. One fluid represents the protons of 
the interstellar plasma component as well as the 
solar wind plasma. The remaining fluids model 
three thermodynamically distinct populations of 
neutral H. Each of the neutral fluids interacts with 
the plasma through resonant charge exchange, us- 
ing the Fite et al. (1962) cross section, and all 
neutrals are subjected to photoionization which 
depends on the squared distance to the Sun. Ra- 



7 



diation pressure is assumed to balance gravity. For 
a detailed description of the numerical model, and 
the underlying physics, see Zank et al. (1996b) and 
Zank (1999). 

Modeling of heliospheric neutrals in the multi- 
fluid code as a superposition of three indepen- 
dent neutral fluids is an approximation to the 
general, non-Maxwellian neutral distribution func- 
tion. Neutrals can also be calculated on a ki- 
netic level without such an approximation (e.g., 
Baranov k Malama 1993; Miiller et al. 2000). Re- 
cent comparison studies (Alexashov & Izmodenov 
2005; Heerikhuisen et al. 2006) have shown that 
the multi-fluid and kinetic methods are in essen- 
tial agreement in most of the flne details concern- 
ing heliospheric geometry and shock locations, and 
even the distributions look similar. For this rea- 
son, we adhere to the computationally less costly 
multi-fluid method for this study. 

For the numerical models, the solar wind at 1 
AU is assumed to be independent of longitude, lat- 
itude, and time, with values of 5.0 cm~^ for the 
plasma density, a temperature of 10^ K, and a ra- 
dial velocity of 400 km s^^. There are two models 
(23 and 24) with slight variations on these solar 
wind conditions (2 x 10^ K and 500 km s"^, re- 
spectively). The models of the interaction of the 
solar wind with the ISM are carried out in a he- 
liocentric frame of reference, and are effectively 
two-dimensional as we assume azimuthal symme- 
try about the stagnation axis (the axis parallel to 
the LISM flow that contains the Sun). To sat- 
isfy this assumption we also neglect heliospheric 
and interstellar magnetic fields. The interstellar 
medium is prescribed as four boundary conditions 
at a suitably large distance from the Sun (a typ- 
ical value is 1000 AU for heliosphercs where the 
interstellar bow shock is less than 500 AU from 
the Sun) . The boundary parameters are the LISM 
H° and H"*" number densities, and the (common) 
hydrogen velocity and temperature. 

3.1. Contemporary LISM 

Inferred values of the contemporary interstellar 
boundary parameters are v = 26.3 km s~^ and 
T ~ 6300 ±340 K (Witte et al. 1996; Witte 2004). 
The contemporary interstellar proton and neutral 
H densities are not well constrained but should 
lie in the range from 0.04-0.14 cm~^ and 0.14- 
0.24 cm~^, respectively (e.g. Zank 1999; Slavin & 



Frisch 2002), and models 4 10 (Table 1) fit within 
these contemporary constraints. We choose model 
5 (n(H+) = 0.047 cm'^, n(H°) = 0.216 cm'^, v = 
26 km s-^, and T = 7000 K) as the highlighted ex- 
ample representing the contemporary heliosphere. 
This model has been described in some detail pre- 
viously by Miiller & Zank (2004), whose Figure 1 
displays two-dimensional maps of plasma temper- 
ature and neutral density, featuring this model's 
heliospheric boundaries together with the neutral 
hydrogen wall. Here, we do not repeat this figure, 
but rather display the plasma temperature (top 
panel) and density (bottom panel) along the stag- 
nation axis as solid lines in Figure 2, together with 
the neutral H density (dash-dotted line). 

The heliospheric boundaries are clearly visible 
as discontinuities in Figure 2, and labels are pro- 
vided next to the temperature profile. The solar 
wind is supersonic at 1 AU, and expands radially 
before undergoing a transition (temperature and 
density increase) at the termination shock (TS). 
The TS is asymmetric, with a nose distance of 
99 AU and a tail distance of 216 AU. In the he- 
liosheath, the region of the shocked and heated so- 
lar wind, the solar wind plasma gets directed tail- 
ward, and is separated from the interstellar plasma 
by the heliopause (HP), with discontinuous den- 
sity and temperature. The stagnation point (the 
nose of the hehopause) is at 148 AU. The LISM is 
supersonic, and consequently there is an interstel- 
lar bow shock (BS) at 285 AU upwind. 

The thermodynamically distinct plasma regions 
define the characteristics of the three neutral flu- 
ids used in the four-fluid model. The compo- 
nent 1 neutral population consists of neutrals from 
the ISM. Tlicy typically are warm and of mod- 
erate bulk speed. Neutrals born through charge 
exchange in the hot hclioshcath between the TS 
and the HP form a different neutral population 
(component 2) which is hot, with correspondingly 
high thermal speeds. The third component con- 
sists of neutrals born inside the TS in the super- 
sonic solar wind; correspondingly, component 3 
neutrals are fast and warm. While the approx- 
imation of the real neutral distribution function 
throughout the heliosphere by a superposition of 
three Maxwellians (the three neutral fluids) intro- 
duces inaccuracies into the models, such a hydro- 
dynamic treatment represents the overall neutral 
H distribution well, as evidenced for example by 
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the successful matching of modeled and observed 
Ly-a absorption by heliospheric H° (Wood et al. 
2000b). 

Figure 2 also shows the temperature and den- 
sity profile of a plasma-only model (model 1, 
dashed lines, TS at 90 AU, HP at 300 AU), which 
will be described in detail in the next section. The 
contrast between the two plasma temperature pro- 
files is due to the effect of charge exchange: The 
pickup process (here, the charge exchange of a 
solar wind proton with an interstellar neutral H 
atom in the supersonic solar wind region) deposits 
energy into the supersonic solar wind and reverses 
the effect of adiabatic cooling (and also slows the 
solar wind), such that the effective solar wind tem- 
perature of model 5 turns upward, whereas model 
1 follows an adiabatic cooling law. In both in- 
ner and outer heliosheath, heat transport by neu- 
tral hydrogen with subsequent secondary charge 
exchange introduces gradual temperature gradi- 
ents between the discontinuities, which are absent 
in model 1 which has no neutrals, no charge ex- 
change, and no anomalous heat transport across 
the HP. 

In the contemporary heliosphere (model 5), 
there is a hydrogen wall between the BS and the 
HP, with a peak density of 0.502 cm^^ = 2.3 n(H°) 
(Figure 2 bottom, dash dot line). The extra wall 
material consists of slower neutral hydrogen born 
from charge exchange with the subsonic interstel- 
lar plasma downstream of the BS (the outer he- 
liosheath). The outer heliosheath plasma is inter- 
stellar plasma that is slowed and heated at the 
BS, and further affected by additional momentum 
loss and energy gain through charge exchange of 
component 2 neutrals that cross from the inner 
into the outer heliosheath. The hydrogen wall is 
accompanied by an elevated plasma density (Fig- 
ure 2 bottom, solid line) because of the plasma 
slowdown. 

The neutral atom density at the upwind TS is 
0.098 cm-3 = 0.46 n(H°), the latter ratio being 
called the filtration factor because it links the in- 
terstellar density to the neutral density in the in- 
ner heliosphere, accounting for the loss processes 
along its path. The neutral density at 5 AU on 
the upwind stagnation axis is 0.036 cm~^ = 0.17 
n(H°). Closer to the Sun, photoionization and 
solar wind charge exchange deplete neutral H ex- 
ponentially. In the tail direction, the stagnation 



axis re-populates slowly with off-axis neutral H. 

Models 4-10 of Table 1 loosely fit within the 
constraints of contemporary observations, and 
generally represent ionization and density levels 
appropriate for low column density ISM where 
both radiation field and abundances of cooling 
trace elements might vary. Table 2 gives the cor- 
responding results for all models, and it can be 
seen that the results from models 4-10 only vary 
modestly. The size of the heliosphere changes, 
with HP locations ranging from 100-150 AU, and 
neutral hydrogen filtration varies from 0.3-0.5, 
while the relative peak wall density remains es- 
sentially unchanged. 

3.2. Hot Local Bubble 

The interior of the Local Bubble void is hot 
and nearly completely ionized. The absence of in- 
terstellar neutral hydrogen in model 1 simplifies 
the heliospheric physics considerably, because in 
all other models the atom-ion process of charge 
exchange generates distinct features in the atom 
distribution function, and also alters the plasma 
by coupling ions and atoms. 

We adopt the ISM parameters of n(H+) = 0.005 
cm-^, n(H°) =0,v = 13.4 km s-\ and log T(K) = 
6.1 (model 1 in Table 1), with the velocity based 
on the Dehnen & Binney (1998) solar apex mo- 
tion since the plasma is assumed at rest in the 
LSR. The speed of sound in such a plasma is 190 
km s~^, and the Sun therefore moves subsonically 
through this medium (Mach 0.07). In this case, 
the isotropic thermal interstellar pressure domi- 
nates the ram pressure, and the termination shock 
is spherical at a distance of 90 AU from the Sun. 
This distance is comparable to that of the contem- 
porary heliosphere. The distance to the nose of 
the heliopause is 300 AU, which makes this sheath 
very large in comparison to the contemporary he- 
liosphere above. The temperature in the sheath 
reaches values as high as 2.2 x 10® K. Figure 2 con- 
tains the plasma temperature profile (top) along 
the stagnation axis, and the plasma density (bot- 
tom), as dashed lines. In the termination shock 
transition, the density jumps by a factor of 3.8, 
and the wind speed decreases to 100 km s~^. 

In this model there are no neutral atoms in 
the entire system of solar wind and interstellar 
medium, and hence the particle content is dif- 
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ferent from that observed today. There are no 
pickup ions (PUI) produced by charge exchange, 
and therefore there are no anomalous cosmic rays, 
and no sfowdown or heating in the supersonic so- 
lar wind beyond the inner solar system. However, 
thanks to its large ram pressure the solar wind 
plasma still provides an effective shield of the so- 
lar system against the million-degree plasma of the 
Local Bubble. 

As a plasma-only model with subsonic inter- 
stellar boundary conditions, model 1 is ideally 
suited for comparison to previous, analytical stud- 
ies of the heliosphere under the assumption of in- 
compressibility of the flow outside of the TS. For 
the simple case of a flow around a rigid sphere 
(as a stand-in for the heliopause), the interstel- 
lar velocity on the stagnation axis behaves as Voo 
(1 — r^p/r'^), where thp is the distance to the 
HP nose, and the uncontaminated interstellar 
flow velocity. This analytical behavior matches 
the one of model 1 very well (Figure not shown 
here). In a detailed treatment, Suess & Nerney 
(1990) calculate the flow streamlines from a re- 
alistic, pressure-balanced TS outward. Their TS 
distance formula (equation (7) below) predicts it 
to be at 81.8 AU, which compares well with the 90 
AU found in model 1, given that the latter (numer- 
ical) model does not presuppose incompressibility. 

It has to be noted that the models do not ac- 
count for the interplanetary magnetic field, as do- 
ing so requires a three-dimensional treatment of 
the problem that is outside the scope of this paper. 
In particular for model 1 in which the mitigating 
aspects of the neutral-plasma interaction are ab- 
sent, the pile-up of magnetic field at the nose of 
the HP (Axford-Cranfill effect; e.g. Nerney et al. 
1993) may be a contributor to the overall pres- 
sure balance in the upwind directions of the inner 
heliosheath. The plasma would not decelerate as 
quickly as l/r^ downstream of the TS, as is ap- 
proximately the case in model 1. Consequently, 
the heliopause would be expected to shift further 
away than the location identified in this paper. 
The unknown interstellar magnetic field strength 
in the Local Bubble has the potential to shift the 
pressure balance of the heliosphere as well. 

3.3. Dense ISM 

Three models in Table 1 represent the effect of 
an ISM that is denser than that of the contempo- 
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Fig. 2. — One-dimensional profiles along the stag- 
nation axis, with the Sun at center and the LISM 
coming from right. Top: plasma temperature 
of model 1 (the Local Bubble case; dashed) and 
model 5 (representing contemporary conditions; 
solid). The hcliospheric boundaries of model 5 
are marked in the plot. The bottom panel con- 
tains the corresponding densities (plasma model 
1, dashed; plasma model 5, solid; model 5 neutral 
H with dash-dot pattern). 




Fig. 3. — Two-dimensional maps of plasma tem- 
perature (left) and neutral H density (right) for 
the high density case 17, with the Sun at center 
and the LISM coming from right. The transition 
from white to medium gray in the plasma temper- 
ature is the interstellar bow shock; the dark shades 
are the hot heliosheath and heliotail. The orbits 
of Saturn, Uranus, and Neptune are sketched as 
dotted hues. 
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Fig. 4. — Similar to Figure 2, 1-D profiles 
of plasma temperatures (top), number densities 
(middle), and parallel velocities (bottom), for 
model 17 (high density case; long and short 
dashed), and model 25 (high velocity case; solid 
and dash-dot). 




Fig. 5. — Two-dimensional maps of plasma tem- 
perature (left) and neutral H density (right) for 
the high speed case 25. Note the triple point at 
about (-105, 30). The orbits of Saturn, Uranus, 
and Neptune are sketched as dotted lines. In the 
neutral density (right panel), the hydrogen wall is 
clearly visible. 



rary heliosphere by a factor of ^50 (models 15- 
17), but not as dense as models investigated by 
Yeghikyan & Fahr (2003) with ntot = 100 cm'^. 
Additionally, models 14 and 27 are denser by a 
factor of 3-4 above the contemporary value. The 
high density leads to a high interstellar ram pres- 
sure, and therefore the resulting heliosphere tends 
to be smaller. 

As an example, Figure 3 displays 2-D maps of 
the hydrogen density and plasma temperature of 
model 17 {n(H°)= 15 cm~'^), and Figure 4 shows 
stagnation axis temperature, density, and velocity 
profiles (dashed). The TS is asymmetric with a 
nose distance of 9.8 AU and a tail distance of 23 
AU (Table 2). The HP is located at 16 AU, and 
the BS at 34 AU upwind. The TS is weak, with an 
upwind compression ratio of 1.8. Atypically, the 
HP is not a sharp temperature gradient as in most 
other models, but the temperature profile is more 
washed out by frequent charge exchange. The bow 
shock has a moderate compression ratio of 3. 

In spite of the small heliosphere, neutrals and 
protons are coupled tightly because of the high 
neutral density, so that the hydrogen wall starts 
immediately downstream of the BS. It has a peak 
density of 3.1 n(H°), but charge exchange is fre- 
quent enough that the neutral density at the 
TS (the filtration factor) is 0.12 n{B°), which is 
among the lowest of the 27 models considered here. 
Even so, the absolute neutral density is quite high, 
leading to a pronounced solar wind slowdown. On 
the upwind stagnation axis, the solar wind speed 
decreases to 260 km s~^ upstream of the TS. In the 
tail region, the solar wind plasma slows to 160 km 
s~^, but frequent charge exchange decreases that 
value in the hcliotail, to a common plasma/neutral 
speed of ^60 km s^^ by ^100 AU downwind of the 
Sun. 

3.4. High Velocity ISM 

Several models in Table 1 test the response of 
the heliosphere on which a high velocity ISM im- 
pinges (models 14 and 18-27). The high velocity 
generates a large ram pressure, making the result- 
ing heliosphere smaller, similar to, but more elon- 
gated than, the high density cases. We present 
model 25 with v ~ 86 km s^^ in more detail here. 
If the corresponding cloud had a thickness com- 
parable to the high velocity ISM towards 23 Ori 
(<0.12 pc, §2.2), it would pass over the Sun in less 
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Table 2 
Model Results 
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233 
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2.3 
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0.23 
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99 


148 
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216 
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0.46 
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0.29 


0.10 


10 


79 


104 


230 


166 


2.3 


0.40 


0.17 


11 


253 


358 


812 


463 


2.1 


0.28 


0.08 


12 


144 


197 


365 


331 


2.4 


0.50 


0.18 


13 


86 


119 


242 


180 


2.3 


0.39 


0.17 


14 


11 


14 


21 


52 


3.3 


0.82 


0.60 


15 


14 


26 


100 


37 


3.0 


0.12 


0.08 


16 


8.2 


12 


23 


28 


7.0 


0.14 


0.10 


17 


9.8 


16 


34 


23 


3.1 


0.12 


0.09 


18 


91 


126 


227 


227 


2.9 


0.35 


0.12 


19 


62 


85 


136 


173 


3.2 


0.50 


0.21 


20 


52 


69 


108 


159 


3.3 


0.59 


0.26 


21 


45 


62 


92 


152 


3.5 


0.68 


0.31 


22 


32 


44 


63 


122 


2.8 


0.97 


0.51 


23 


32 


44 


63 


121 


2.9 


0.96 


0.52 


24 


39 


54 


79 


139 


3.0 


0.82 


0.42 


25 


26 


34 


50 


112 


2.7 


1.00 


0.64 


26 


38 


44 


77 


141 


5.3 


0.87 


0.25 


27 


21 


31 


46 


72 


3.4 


0.88 


0.54 
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than ^1400 yr. Figure 5 displays 2-D maps of hy- 
drogen density and plasma temperature of model 
25, showing the heliospheric boundaries and fea- 
tures clearly. The 1-D stagnation axis profiles in 
temperature, density, and velocity are displayed in 
Figure 4. 

In most models, the TS, taken as a 3-D smface, 
is spherical (model 1), or nearly spherical with 
an upwind/downwind asymmetry. These cases 
are characterized by a heliosheatli and heliotail 
plasma that are subsonic throughout. In contrast 
to this, the shape of the TS of a high velocity 
ISM heliosphere, such as model 25, is qualitatively 
different, now resembling a rocket shape. The 
initially subsonic plasma at the nose of the he- 
liosheath accelerates to supersonic speeds in the 
nozzle-shaped region between the TS and the HP. 
However, to match the subsonic heliotail plasma 
and the supersonic hclioshcath plasma requires 
both a shock to decelerate the flow and a tangen- 
tial discontinuity to adjust the density. A char- 
acteristic triple point occurs where heliosheath 
shock, termination shock, and the tangential dis- 
continuity meet. Figure 5 shows this morphology 
in the example of model 25. 

The TS is highly asymmetric, with a nose dis- 
tance of 26 AU and a tail distance of 112 AU 
(Table 2). The upwind TS compression ratio s 
is s = 2.9. The HP is at 34 AU, and the BS at 
50 AU upwind. The bow shock is quite strong, 
with a post shock plasma speed of 24 km s~^ 
(Figure 4, bottom panel), temperature of 10^ K 
(top), and a compression ratio of 3.4 (middle). Be- 
cause of the large neutral velocity in the post-bow 
shock rc^gion, the neutral mean free path (mfp) for 
charge exchange is initially ^30 AU, larger than 
the outer heliosheath, and shortens only gradually 
as the effective neutral velocity decreases to 31 km 
s~^. Consequently, the model 25 hydrogen wall 
between the BS and the HP is not very thick, but 
reaches a peak density of 2.7 n(H°) about ~11 AU 
downstream of the BS. This distance from the BS 
is of the same magnitude as in the more moderate 
cases, however, in the high speed case, it brings 
the peak close to the HP already. As the TS is so 
close to the hydrogen wall, the filtration factor is 
1.0, i.e. the neutral density at the TS equals that 
of the LISM (and is still 0.6 n(H°) at 5 AU, Figure 
4, middle panel). These filtration factors close to 
unity (models 14, 22 - 27 in Table 2) seem only 



possible when a high interstellar velocity combines 
with a modest or low density so that the peak hy- 
drogen wall occurs close to the HP without room 
for depletion of neutral H between peak and HP. 
In the similar sized dense heliosphere (§3.3), the 
charge exchange mfps are shorter, the peak hydro- 
gen wall is attained farther away from the HP, and 
charge exchange upwind close to the HP spreads 
the H flow and leads to a density decrease already 
before the H° flow crosses the HP (Figure 4). 

Finally, we note that for the high speed mod- 
els 14 and 22 - 27, the original numerical grid 
with a 5° angular resolution is too coarse, lead- 
ing to errors during the transport along directions 
with high velocity components. In these cases, we 
choose a grid with a 2° angular resolution, which 
cures the problem so that pre-BS density values 
on the stagnation axis are within 10% or better of 
the LISM values. They are always more accurate 
in off- axis directions. 

4. MODEL CORRELATIONS 

In addition to the individual model results dis- 
cussed in the preceding section. Table 2 contains 
the same key results characterizing the boundary 
locations and neutral H content for all 27 mod- 
els. The model numbers refer to the correspond- 
ing boundary parameters listed in Table 1. Model 
1 stands out from the rest in that neutrals are ab- 
sent. For the subsonic models 1 and 2 the LISM 
pressure is dominated by the thermal pressure. All 
other models are ram-pressure dominated. This is 
among the reasons why the results discussed below 
cover only a subset of the vast parameter space. 

4.1. Plasma structure 

When comparing the 27 models, an obvious re- 
sult is the variation in the size of the heliosphere, 
as expressed in the location of upwind TS, HP, and 
BS in Table 2, as well as the distance of the down- 
wind termination shock (TS^). These distances 
are set by balancing the solar wind pressure and 
interstellar pressure (e.g. Holzer 1989). A small 
heliosphere is caused by a large LISM pressure. 
Models 2, 3, 11 and 12 are especially large due to 
a lower LISM plasma ram pressure, namely, a low 
LISM velocity in the case of models 2, 3 and 11, 
and a low density for models 11 and 12. Over- 
all, the upstream distances from the Sun to the 
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TS range from 8 to 260 AU, and those to the HP 
range from 12 to 400 AU (see Table 2). The bow 
shock of a cool, slow, tenuous LISM (model 11) is 
as far as 810 AU away from the Sun. 

Because the overall system is pressure balanced, 
the locations (heliocentric distances) of TS, HP, 
and BS arc correlated with one another. Taking 
all models except model 1, and adding the results 
of another systematic (as yet incomplete) param- 
eter study with v — 26.24 km s"-"-, the correlation 
between rxs, the distance of the upwind TS, with 
the distance of the upwind heliopause thp , is 

rnp = (1.39 ±0.01) TTS, (1) 

obtained with a linear regression analysis after as- 
cribing uncertainties to rxs and rnp due to grid 
resolution and HP stability. The intercept value of 
the analysis is consistent with zero. The omitted 
model 1 ratio rHp/rxs = 3.3; the other subsonic 
case 2 is an outlier with a ratio of 1.55. The filled 
circles in Figure 6 represent each model as a (TS, 
HP) pair, and the corresponding straight line is 
the linear fit of equation (1). Model 1 is addition- 
ally marked with an " x " . 

For models with a supersonic LISM, the upwind 
bow shock TBS is located about twice as far as the 
heliopause, 

TBS = (1-95 ±0.05) THP. (2) 

A direct correlation between TS and BS yields 
''BS = (2.70±0.09) rxs- In deriving the BS correla- 
tions, models 3 and 11 were excluded as the largest 
models; the above relation underprcdicts the BS 
distance of model 3 and 11 by 81 AU and 114 
AU, respectively. Figure 6 shows the model (TS, 
BS) pairs as triangles. In contrast to the excellent 
HP data fit, the BS locations are more scattered 
around their straight line fit. For the smaller he- 
liospheres, the BS lies systematically more inward 
than predicted from the TS location via the pro- 
portionality fit (2), and apparently obeys a differ- 
ent linear relation. 

For the interstellar ram-pressure dominated 
models considered here, the upwind- downwind 
asymmetry of the TS is essentially constant. Ex- 
cluding the rocket-shape models that have a triple 
point in their TS, the relation for the downwind 
distance rxsd of the TS is 

rxsd = (2.08 ±0.04) rxs. (3) 



The appearance of a triple point in the remain- 
ing models clearly yields a different dependence, 
namely, 

rxsd = (1.35 ± 0.12) rxs ± (82 ± 7)AU. (4) 

Again, these fits are shown in Figure 6, with open 
circles representing their corresponding model 
data. Both subsonic models 1 and 2 are out- 
liers and were omitted. They obey a different 
asymmetry law than the rest of the models, stem- 
ming from the qualitatively different pressure dis- 
tribution along the HP in the absence of a bow 
shock. Model 11 is an outlier with an rxsd smaller 
by 63 AU than predicted by the above relation. 
For small heliospheres, relation (3) fits the high 
density cases 15-17 well, but underprcdicts the 
asymmetry of models 14 and 27, which due to their 
high velocity might be on their way to the rocket- 
shaped cases. The mentioned outliers are not 
included in the factors of equations (3) and (4). 
The pressure distribution in the rocket-shaped 
cases goes hand-in-hand with a more inward loca- 
tion of the bow shock; the models where the BS 
location is inconsistent with Eq. (2) are the ones 
where the heliosphere is of that shape. 

In the supersonic solar wind, the total pressure 
is dominated by the plasma ram pressure. Since 
most of the kinetic energy of the solar wind is 
converted into heat at the termination shock, it 
is effectively the upstream solar wind ram pres- 
sure that balances the total interstellar pressure 
-Plot = -Ppi + Ppv'^ + + Pnv^, the sum of plasma 
thermal and ram pressure and neutral H thermal 
and ram pressure. Hence a simple 1-D pressure 
balance would be achieved, assuming constant so- 
lar wind velocity vsw and an dependence of 
the density for the supersonic solar wind region, 
at the radial distance of 

Tpb/ri = VPi/Ptot (5) 

where Pi = PiVg^^ is the solar wind ram pressure, 
and pi the solar wind density, both taken at ri = 1 
AU. This basic pressure balance distance assumes 
that all kinetic energy is converted into heat at 
the TS. Improving on this by using the Rankine- 
Hugoniot relations together with treating the he- 
liosheath and interstellar flows as incompressible, 
and assuming the ISM to be at rest in the helio- 
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centric frame, the TS is calculated to be at 



7 + 3 7 

= ]l n:^) ^ y 8 

(Zank 1999), where 7 is the ratio of specific heats, 
and is set to 5/3 in the second equality as well as 
for all the models in this paper. 

Suess & Nerney (1990) have calculated the case 
of an ISM that is moving with respect to the Sun, 
in which the nose TS position calculated from 
equation (6) is corrected by a weak dependence 
on the ratio of interstellar (v) and solar wind ve- 
locities, 

rTS2 = \ — — r rpb — 1 — 

16 1 f,_J^Y (T 

~ 5 V vl^J ■ 

For 7 = 5/3, the numerical factors of equations (6) 
and (7) agree with each other to within 0.4%, and 
the extra velocity-related factor of (7) represents a 
non- negligible reduction in rTS2 compared to ttsi 
only for high velocity ISM cases, such as model 14 
of this paper (a 8% reduction in this case). 

A number of assumptions enter into the deriva- 
tion of equations (6) and (7). Among them is that 
these formulas are only valid for a subsonic ISM, 
with neutral H absent, and Ptot should be domi- 
nated by the ISM thermal plasma pressure. How- 
ever, the validity of (6) and (7) can be extended 
to the supersonic LISM case, because an interstel- 
lar bow shock will convert the flow to subsonic 
speeds. It is therefore possible to admit the sum 
of all plasma pressure contributions into Ptot of 
equation (5). If neutral H were tightly coupled 
to the plasma (i.e. if the mean free paths were 
very short compared to typical heliospheric length 
scales), then the neutral H pressure contributions 
can justifiably be included in Ptot as well. How- 
ever, the neutral-plasma coupling is neither zero 
nor very strong, so that the solar wind/LISM pres- 
sure balance for the heliospheres modeled in this 
paper is more complicated. 

The 27 models show a correlation between rTS2 
and the upwind TS location (and consequently, 
the upwind HP and BS, equations 1 and 2). We 
present here two of the rTS2 distance estimates. 



One {vpi) is obtained by using solely the LISM 
plasma pressure in equations (5, 7) (Ptot, pi = 
Ppi + PpV^), effectively neglecting neutral H alto- 
gether. The altc;rnativc distance prediction rmin 
uses Ptot = Ptot.pi + Ph + PHf ^, taking the entire 
pressure contribution of neutral H into account as 
well. The correlations are 

TTS = (0.71 ± 0.08) rpi; Tp, = rTS2(Ptot,pi)(8) 
tts = (1.35 ± 0.02) r„ji„; rmi„ = rTS2 (Ptoi(P) 

The data and the regression fits arc shown in Fig- 
ure 6, top, with open diamonds and squares, re- 
spectively. The plasma-only correlation (8, dia- 
monds) is poor, and especially pairing models 3 
with 11, and 26 with 27, shows that the same 
plasma total pressure does lead to different he- 
liospheres depending on the neutral contribution. 
In addition, the high density cases 15-17 are com- 
plete outliers (as are 26 and 27) where the plasma 
pressure alone overpredicts the heliospheric size. 
For the correlation (9) (squares), it is not neces- 
sary to exclude any of the supersonic heliospheres, 
in particular not models 15-17. Only model 2 is 
an outlier, but not model 1. 

The predicted and modeled TS distances in the 
plasma-only relation (8) scatter aroimd the fit too 
much, and neutral dominated models are excluded 
outright. The deviations from the predicted loca- 
tions indicate changes to the pressure balance due 
to charge exchange with neutral H. The pressure 
balance gets shifted inwards by pickup ion produc- 
tion in the supersonic solar wind that reduces the 
supersonic wind speed and hence its ram pressure. 
On the other hand, the TS pressure balance can 
get shifted outwards by other effects, among them 
the deceleration of the outer heliosheath plasma 
by charge exchange with secondary neutrals. Re- 
lation (9) has a much smaller scatter and therefore 
is better suited to predict rxs , with the caveat that 
the neutral-plasma interaction drives the neutrals 
out of equilibrium, and pressure balance can only 
be described by an empirical factor of 1.35. 

This study is focusing solely on the variation 
of interstellar boundary parameters and the helio- 
spheric response. While we do not wish to study 
the question of different solar wind conditions on 
the heliosphere here, it is interesting to note that 
all the results so far only depend on Pi, the solar 
wind ram pressure at 1 AU, and that the loca- 
tions of TS, HP, and BS are dictated by a pres- 
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Fig. 6. — Correlation between the upwind ter- 
mination shock location tts a-nd other distances, 
including, in the bottom panel, the upwind he- 
liopause HP (red circles), the upwind bow shock 
(green triangles), and the tailward TS distance 
(open blue circles), which are all model results. 
Model 1 results are marked with "x". The up- 
per panel shows the TS correlations with the dis- 
tance rTS2 theoretically derived from plasma pres- 
sure balance (pink diamonds) and from total pres- 
sure balance (blue squares), together with empiri- 
cal linear fits. Note the aspect ratio of 2:1 in both 
panels. 



sure balance between solar wind and LISM, with 
empirical correction factors. In reality, the solar 
wind is time dependent on an 11 year cycle and on 
smaller, episodic timescales. The variation in so- 
lar wind ram pressure leads to small variations in 
the heliospheric boundary locations (e.g., Zank & 
Miiller 2003; Izmodenov et al. 2005) that are qual- 
itatively consistent with Eq. (9) . Similarly, the in- 
creased ram pressure in polar directions from the 
fast polar wind during solar minimum results in 
larger heliospheric distances in these directions as 
compared to isotropic slow solar wind (Pauls & 
Zank 1997; Tanaka & Washimi 1999). An exam- 
ple of the reaction of the heliosphere to different 
solar wind ram pressures while holding the ISM 
environment constant is the comparison of mod- 
els 22-24 in Table 2, where the boundaries move 
outward for a ram pressure increase (model 24), 
but do not change for a doubling of the solar wind 
thermal pressure (model 23). 

Recent investigations of astrospheres around 
other cool main-sequence stars (Miiller et al. 2001; 
Wood et al. 2002, 2005a) use a range of stellar 
wind ram pressures Pi that are different from the 
solar wind, but a modeling strategy identical to 
the one for the heliosphere. A cursory analysis of 
these models (plots not shown here) shows that 
equations (1) - (9) hold also for all these cases 
within the stated accuracy. This result under- 
scores the argument of pressure balance, and ex- 
tends the results of this section to other ram pres- 
sure regimes such as for astrospheres carved out 
of the ISM by coronal stellar winds. 

4.2. Heliospheric neutral hydrogen 

The density 71^5(11) of neutral hydrogen that 
crosses the termination shock at the upwind stag- 
nation axis varies over a considerable range in the 
27 models, from 0.01 to 0.33 cm~^, with four larger 
values between 0.9 and 2 cm~^ for the high den- 
sity models 15-17 and 27. The filtration ratio / 
is the neutral density at the TS divided by the 
interstellar neutral density n(H°), well upstream 
of the BS (to avoid contamination by component 
2 and 3 neutrals), and these relative values vary 
from 0.1 to 0.6, with a few higher values 0.7-1.0 
for higher velocities. The filtration is listed as /ts 
in Table 2 for all 27 models, as are similar ratios 
for the peak hydrogen wall density (/peak) and for 
the density at 5 AU on the upwind stagnation axis 
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Fig. 7. — Correlations between the neutral H den- 
sity at 5 AU upwind, and at the TS (circles) 
and the peak density inside the hydrogen wall, 
max{nH (wall) ) (diamonds) . 
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Fig. 8. — Correlation between the filtration ratio 
/=nTs(H)/n(H'^) and the interstellar velocity v 
(circles), and correlation of the 5 AU density ra- 
tio n5Au(H)/n(H°), with v (diamonds). The solid 
lines are linear fits while the plot itself is double- 
logarithmic. 



(/sAu)- We choose 5 AU as a fixed reference dis- 
tance with the expectation that photoionization is 
not yet important at this distance. 

The neutral hydrogen at the TS is comprised 
of original interstellar material, and of slower sec- 
ondary neutrals created upwind of the HP which 
form the hydrogen wall. After crossing the HP 
these neutrals get depleted in the heliosheath by 
charge exchange which replaces them with neu- 
trals mostly in outward directions. For high veloc- 
ity models, this loss process starts from a higher 
HP density so that higher filtration ratios occur 
(§3.4). Both the filtration ratios and the absolute 
TS neutral densities correlate well with the neutral 
density n5Au(H) at 5 AU. The correlations are, 

nTs(H) - (1.43 ± 0.02) n5Au(H)°-''" °' (10) 

/ - (1.55 ±0.13) [n5Au(H)/n(H")]°-''t'i^?1 

where all models have been included in (10), and 
the high density models 15-17 had to be excluded 
as outliers in (11) (plot not shown). The data 
and the fit (10) are displayed in Figure 7 by circle 
symbols. This correlation is interesting in that it 
relates a density at a fixed distance to a density at 
the TS regardless of the actual distance of the TS. 
The net loss in the region of the supersonic solar 
wind is not sensitive to the length of the neutral 
particle trajectories between the TS and 5 AU be- 
cause most of the charge exchange relevant to the 
5 AU density takes place immediately upwind of 
that distance, given that the local plasma density 
is larger with smaller distance, resulting in smaller 
charge exchange mean free path lengths. 

The peak densities in the hydrogen wall range 
from 1.1 n(H°) to massive hydrogen walls of 7 
n(H°), corresponding 0.09-105 cm~^ in absolute 
units. The absolute peak density is weakly corre- 
lated to the 5 AU value, as demonstrated by the 
diamond symbols in Figure 7. Even on a log-log 
plot, however, a large scatter of the results around 
possible power laws is evident. According to a re- 
gression analysis that omits the high density mod- 
els 15-17, 



,(H) = (7.1 ±1.5) n5Au(H) 



0.81±0.06 



(12) 



Nonetheless, the fit is poor. The same problem is 
encountered when relating the neutral results to 
the interstellar velocity. There is a general trend 
of higher neutral densities with higher velocity. 
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Figure 8 shows the filtration results (circles) and 
the relative densities at 5 AU (diamonds) plot- 
ted against the interstellar velocity. It is evident 
that any correlation derived from these model re- 
sults is not unique, as there are many data points 
around ?; « 26 km whose normalized neu- 
tral results vary by a factor of 2.7. The fit lines 
shown in Figure 8 are / = v/{73 ± 7 km s~^) and 
n5Au(H)/n(H°) = v/{U0±8 km s'^). 

By way of a cautionary note, it should be noted 
that the above relations are likely model depen- 
dent, in that using another self-consistent model- 
ing strategy for the heliosphere, such as a parti- 
cle kinetic model for the neutral H (Baranov & 
Malama 1993; Miiller et al. 2000; Heerikhuisen 
et al. 2006), might result in different coefficients 
for the relations (l)-(4), and (8)-(12). Short of 
carrying out a parameter study with these alter- 
nate models and comparing the results, it is im- 
possible to incorporate this systematic error into 
the error estimates of the coefficients given in the 
above relations. 

5. DISCUSSION 

5.1. Heliospheric Morphology and Neu- 
trals 

We are now in a position to discuss some of 
the consequences that different interstellar envi- 
ronments have for the solar system and for Earth. 
There are several models listed in Table 2 rep- 
resenting heliospheres that are so small that the 
outer planets find themselves beyond the super- 
sonic region of the solar wind, at least for parts of 
their orbits. As examples. Figures 3 and 5 have 
the orbits for the gas giants Saturn, Uranus, and 
Neptune marked in them. An assumption in these 
plots is that the ISM flow vector is close to the 
plane of the ecliptic, as is the case in the current 
environment. For larger angles, the ecliptic TS 
distances are farther due to the upwind/downwind 
asymmetry of the heliosphere. 

In the high density case 17 (Figure 3), Uranus 
periodically crosses the TS, traverses the hot inner 
heliosheath, crosses the HP, and spends part of its 
orbit in the shocked LISM (the outer heliosheath), 
before crossing the HP again in a reverse course 
of events. Neptune is always surrounded by hot, 
shocked plasma (either the solar wind heliosheath 
and heliotail or the interstellar plasma of the outer 



heliosheath), and is never upstream of the TS. 
Similarly, in the high velocity case of model 25 
(Figure 5), part of Neptune's orbit is in the sub- 
sonic solar wind of the heliosheath, and Uranus 
periodically comes close to the TS. 

For planets in the hot inner heliosheath it can 
be expected that solar wind injection into plan- 
etary magnetospheres is more efficient due to in- 
creased thermal plasma velocity and solar wind 
turbidcnce. Also the solar wind magnetic field 
strength is enhanced, and an increased flux in 
both galactic cosmic rays and heliospheric ener- 
getic particles is to be expected. Planets that cross 
into the outer heliosheath will be exposed to not 
quite as hot a plasma, carrying a basically unde- 
termined magnetic field. However, the planet will 
then be exposed to the increased neutral density 
of the hydrogen wall. This neutral flow will strike 
the planetary atmosphere unimpeded and lead to 
atmospheric drag and other effects (Yeghikyan & 
Fahr 2004b). 

None of the 27 models considered here leads to 
such dire predictions for Earth's orbit, and extrap- 
olating from equations (7) and (9) beyond their 
validity, interstellar densities of 1500 cm~^ with 
moderate interstellar velocities, or a velocity of 
345 km s^^ with moderate interstellar densities, 
would be needed to place the TS at 1 AU. It should 
be expected that both these numbers are only a 
qualitative estimate, and that detailed modeling 
of such heliospheres would have to take into ac- 
count additional physical processes (Yeghikyan & 
Fahr 2003, 2004a). 

As many models in this paper have an inner 
boundary at or beyond 1 AU, we want to focus 
on the reference distance of 5 AU instead to as- 
sess the change in particle environment that oc- 
curs with different Galactic environments under 
the assumption of an unchanging solar wind. At 
this distance, the solar wind is supersonic for all 
27 models. The ratio of interstellar (slow) neutral 
hydrogen to solar wind protons at 5 AU upwind 
is between 7% and 25% for the contemporary he- 
liosphere (models 4 10). As already mentioned in 
§4.2, this density ratio scales with interstellar ve- 
locity, ranging from from 2% to 37% (models 2 
and 18-24, respectively; model 14: 120%, model 
25: 45%), and is highest for the high density mod- 
els 14-17 and 27. In particular, there is more than 
seven times more neutral H than solar wind pro- 
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tons at 5 AU in models 16 and 17. 

Naturally, the increased presence of interstellar 
H° increases the rate of charge exchange and hence 
the production rate of fast component 3 neutrals, 
the so-called neutral solar wind (NSW). The NSW 
density at 5 AU for the contemporary heliospheres 
is KiA xlO^'' cm~3, but increases to 0.014, 0.018, 
and 0.023 cm~^ for models 15, 16, and 17, respec- 
tively. A similar drastic relative NSW increase 
should occur at Earth orbit as well for these high 
density cases. In contrast, the high velocity case 
25 only yields a five-fold increase of NSW to 0.002 
cm~^ at 5 AU. 

5.2. Cosmic ray transport model 

As an application of the results discussed above 
we make use of a cosmic ray transport model for 
three demonstrative solutions of GCR phase space 
density. They are calculated for the interstellar 
environments corresponding to the Local Bubble, 
the Lie, and a dense cloud of mostly neutral hy- 
drogen, based approximately on models 1, 9, and 
15, respectively. The cosmic-ray transport model 
we use is discussed in detail in Florinski & Zank 
(200G) and the reader is referred to that paper for 
a complete description. Briefly, the plasma flow 
background obtained from the multi-fluid code is 
used to calculate all 3 components of the helio- 
spheric magnetic field in the azimuthal plane from 
Faraday's law combined with the zero divergence 
condition with a Parker spiral field specified at 
the inner boundary (Florinski et al. 2003b). The 
modified field component of Jokipii & Kota (1989) 
that alters the transport parameters at high he- 
liographic latitude is also included in the model. 
The interstellar magnetic field plays no role in this 
model because it fluctuates on scales that are typ- 
ically much larger than diameter of the cosmic ray 
gyroorbit and hence has little effect on the parti- 
cle's trajectories. Next, magnetic turbulent energy 
{6B^) and the associated turbulence correlation 
length in the solar wind region are computed from 
the hydrodynamic model of incompressible tur- 
bulence transport (Zank et al. 1996a; Matthaeus 
et al. 1999). The model assumes that the num- 
ber of waves propagating parallel and antiparallel 
to the mean magnetic field are equal and ignores 
certain wave propagation effects by neglecting the 
Alfven speed compared with the mean plasma ve- 
locity in the solar wind. 



Little is known observationally about the tur- 
bulent content of the inner heliosheath, and turbu- 
lence transport in that region is poorly understood 
at present. Here we use a simple assumption that 
the turbulent ratio {SB^)/B'^ and the correlation 
length are both constant across the termination 
shock and in the heliosheath. This assumption is 
based on the physics of Alfven wave transmission 
through a quasi-perpendicular shock that yields 
the expression relating the turbulent ratio on the 
two sides of the shock as (McKenzie & Westphal 
1969) 

{SBi) _ {s + l){5Bf) 
Bl ~ 2s Bl ' 

which is not too different from 1 for shocks of mod- 
erate strength (compression ratios s = 2.5-3.0). 

The turbulent content of the solar wind is 
strongly infiuenced by the process of Alfven wave 
generation by pickup ions as they scatter from 
the initial ring-beam distribution onto a bisphere 
in velocity space (Williams & Zank 1994; Isen- 
berg et al. 2003). Because pickup ions are pro- 
duced in charge transfer collisions between so- 
lar wind protons and interstellar hydrogen atoms, 
changes in neutral density are the principal source 
of turbulence variability. The latter translates 
into variations in the amount of GCR modulation 
through an appropriate diffusion model. We use 
the Quasi-Linear theory (QLT) (e.g., Jokipii 1966) 
for the parallel component of the diffusion tensor 
and the Nonlinear Guiding Center (NLGC) the- 
ory (Matthaeus et al. 2003) for the perpendicular 
component. The former is governed by the fiuc- 
tuations with wavevectors oriented parallel to the 
mean magnetic field (the slab component), while 
the latter is determined by the fiuctuations orthog- 
onal to the field (the 2D component). We assume 
that slab fluctuations comprise 10% of the total 
energy with 2D making up the rest, which is in 
agreement with the observed ratio (Bieber et al. 
1996). 

Both QLT and NLGC require a specification 
of the reduced (one-dimensional) turbulent spec- 
trum. The power spectrum measured in the so- 
lar wind at low latitudes consists of a flat en- 
ergy range followed by a Kolmogorov inertial range 
with k~^^^ (e.g., Bieber et al. 1994), while some 
observational evidence points to a dependence 
at high latitudes (Horbury & Balogh 2001). The 
low end of the turbulence spectrum is likely to 



19 



be populated by structures (shocks and disconti- 
nuities) that may be responsible for the modula- 
tion of very high energy (above 1 GeV) particles. 
Accordingly, we use two forms of the power spec- 
trum, identical in the inertial range, but having 
a different spectral index in the energy range: 
(diffusion Model I) and —1 (diffusion model 11). 
As discussed in Florinski & Zank (2006), the first 
model emphasizes modulation in the hclioshcath 
region, while model II is dominated by solar wind 
modulation. Because the amount of modulation 
in the heliosheath is not known at present we use 
two plausible models to cover the full range of pos- 
sibilities. 

5.3. Cosmic ray modulation results 

The different sizes of the heliosphere as well as 
the different particle distributions will affect the 
modulation of GCRs as they pass from the ISM 
through the heliosphere into interplanetary space 
and to Earth. The GCR transport coefficients de- 
pend on the heliospheric magnetic field as well as 
on the level of plasma turbulence. While the mag- 
netic field geometry, its strength, and the proper- 
ties of solar wind turbulence (especially in the he- 
liosheath) are far from having been studied conclu- 
sively, theoretical models of cosmic-ray transport 
can be constructed on the basis of known physics 
(Florinski et al. 2003b; Florinski & Zank 2006). 
Here we perform a series of computer simulations 
of GCR modulation in the global heliosphere for 
models 1 (heliosphere in the Local Bubble), 9 (con- 
temporary heliosphere), and 15 (a high density 
cloud encounter). 

A thicker heliosheath such as in model 1 could, 
in principle, be expected to yield stronger modula- 
tion of GCRs (i.e. a lower particle flux arriving at 
Earth). However, the level of magnetic field turbu- 
lence production, which in turn depends partly on 
PUI production, can be less efficient in low density 
environments. Figure 9 shows the high energy pro- 
ton GCR spectra for the three cases (models 1,9, 
and 15) calculated with diffusion models I and II. 
The right panel of this figure shows particle inten- 
sity at the termination shock thus demonstrating 
the effect of heliosheath modulation. Because low 
energy GCR protons typically do not reach 1 AU, 
we focus on the high energy end of the spectra. 

It follows from our model calculations that the 
GCR environment at Earth for the Local Bubble 
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Fig. 9. — GCR proton differential spectra at 1 AU 
(top) and just beyond the termination shock (bot- 
tom) in the apex direction for the three interstellar 
environments represented by models 1,9, and 15. 
The ISM GCR spectrum assumed at the exter- 
nal boundary (Ip & Axford 1985) is shown with 
a solid black line. The roman numerals following 
the model number refer to turbulence evolution 
models I and II. 1 AU spectra observed during the 
time around the 1995 solar minimum are shown for 
comparison. The low energy data are from IMP-8 
satellite observations (McDonald 1998) while high 
energy data are from BESS balloon measurements 
(Sanuki ct al. 2000). 
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scenario (model 1) would have been less intense 
than at present if diffusion model I was correct, 
or, surprisingly, more intense with model II diffu- 
sion coefficients. The latter is a consequence of the 
large GCR mean free path predicted by the second 
diffusion model for the solar wind region in model 
1 in the absence of PUI turbulence driving. In 
both cases the cumulative heliosheath GCR mod- 
ulation is more important for the Local Bubble en- 
vironment than for the contemporary hcliosphcrc. 
This is a combined effect of a stronger TS, a larger 
decrease in the radial mean free path across the 
shock, and a thicker heliosheath. 

For an encounter with a high density cloud 
(model 15), the heliospheric GCR shielding is 
much less effective than in the contemporary he- 
liosphere, such that the high energy part of the 
GCR spectrum approaches the one assumed for 
the pristine LISM. Here, the opposing effects of 
a relatively large heliosheath diffusion coefficient 
combined with a much smaller extent of the mod- 
ulation cavity, and enhanced PUI turbulence driv- 
ing in the solar wind result in a significant reduc- 
tion in the heliospheric shielding of GCRs. The 
predicted cosmic-ray intensity increase at Earth 
is between 1.4-2.4 (model I) and 4.1-7.6 (model 
II) in the energy interval between 300 MeV and 1 
GeV. 

In addition to GCRs, the distribution of anoma- 
lous cosmic rays which are accelerated at the TS, 
will depend on the background neutral density, 
and the strength and distance of the TS (Florinski 
et al. 2003a). Together, the cosmic ray environ- 
ment at Earth influences the terrestrial magneto- 
sphere as well as climate, atmosphere, and biol- 
ogy (e.g. Scherer 2000; Yeghikyan & Fahr 2004b; 
Prisch et al. 2005). 

6. CONCLUSIONS 

On its path through the galaxy, the Sun has 
encountered (and will encounter) different inter- 
stellar environments. This motivates a parameter 
study to investigate the response of the heliosphere 
to these changing conditions under the assumption 
of a constant solar wind. For conditions that are 
not too far from the contemporary LISM environ- 
ment, the following findings emerge from analyz- 
ing 27 self-consistent multi-fiuid models. 

1. Allowing generous assumptions about the 



Lie morphology, the LIC column density to- 
wards nearby stars indicates the Sun first en- 
countered the LIC gas within the past 40,000 
yr, and the CLIC within the past -60,000// 
yr (where / represents the fraction of space 
filled by the CLIC). The Sun is expected to 
exit the LISM gas cloudlet, which is charac- 
terized by the common LIC velocity, some- 
time within the next ~0-4000 yr. In gen- 
eral, passage through interstellar clouds will 
lead to variations in the heliosphere bound- 
ary conditions over timescales possibly as 
short as lO'' yr. Nearby ISM generally re- 
sembles low column density ISM observed 
elsewhere. 

2. The size of the heliosphere is determined 
by the balance of solar wind and interstel- 
lar pressure. For the investigated parame- 
ter range, in which the LISM is mostly ram- 
pressure dominated, the upwind termination 
shock distance can be estimated by equation 
(9), using equations (5) and (7). This rela- 
tion is derived from a pressure balance ar- 
gument modified by an empirical factor ex- 
pressing the efficiency of the neutral pressure 
contribution to the overall interstellar pres- 
sure. 

3. Heliocentric distances of interest such as the 
heliopause, the bow shock, or the upwind 
and downwind termination shock scale lin- 
early with each other (e.g. thp = 1-39 rxs, 
Tbs = 1-95 Thp)- Therefore, when the up- 
wind termination shock distance is predicted 
in an absolute way as above, the other dis- 
tances can be predicted as well. However, 
the scalability and predictability of the he- 
liosphere size with these relations are only 
applicable to parameter sets in which the 
LISM fiow is supersonic. The subsonic cases, 
when the Sun is surrounded by hot plasma 
or alternately when the Sun and the sur- 
rounding interstellar cloud are comoving in 
space, obey a different set of correlations, 
and are generally more difficult to model nu- 
merically. 

4. For low interstellar velocities, the heliosheath 
and heliotail plasma are subsonic through- 
out, and the ratio of downwind to upwind 
termination shock distance (TS asymmetry) 
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is 2.1. For higher velocities, the heliosphere 
assumes a rocket shape, with a modified 
pressure balance in the downwind directions. 

5. Neutral hydrogen results such as the filtra- 
tion ratio, the peak hydrogen wall density, or 
the density at 5 AU upwind of the Sun, cor- 
relate with each other. Their absolute value 
is weakly correlated to the interstellar veloc- 
ity with which the neutrals arrive at the re- 
spective heliosphere, as the charge exchange 
mean free path depends on this velocity, and 
higher velocities shorten the heliocentric dis- 
tances to the heliospheric boundaries. 

6. For encounters with a high density interstel- 
lar cloud (-^15 cm~^, about 50 times the 
contemporary value), the particle fluxes ar- 
riving at Earth orbit, including interstellar 
neutrals, neutral solar wind, and cosmic rays 
will increase markedly. These changes po- 
tentially affect Earth's atmosphere and its 
climate. The changes in particle fluxes just 
due to a higher interstellar velocity are less 
pronounced. 

7. For the period when the Sun was embed- 
ded in the Local Bubble, particle fluxes were 
reduced substantially. Secondary particles 
like anomalous cosmic rays and neutral solar 
wind were entirely absent, and the galactic 
cosmic ray flux arriving at Earth was com- 
parable to the contemporary flux, or even re- 
duced, depending on the modulation model. 
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